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ABSTRACT 

A Spitzer Space Telescope survey in the NOAO Deep- Wide Field in Bootes 
provides a complete, 8 /im-selected sample of galaxies to a limiting (Vega) mag- 
nitude of 13.5. In the 6.88 deg 2 field sampled, 79% of the 4867 galaxies have 
spectroscopic redshifts, allowing an accurate determination of the local (z < 0.3) 
galaxy luminosity function. Stellar and dust emission can be separated on the 
basis of observed galaxy colors. Dust emission (mostly PAH) accounts for 80% 
of the 8 /im luminosity, stellar photospheres account for 19%, and AGN emis- 
sion accounts for roughly 1%. A sub-sample of the 8 /im-selected galaxies have 
blue, early-type colors, but even most of these have significant PAH emission. 
The luminosity functions for the total 8 /im luminosity and for the dust emission 
alone are both well fit by Schechter functions. For the 8 /im luminosity function, 
the characteristic luminosity is uL*(8.0 /im) = 1.8 x 10 10 L & while for the dust 
emission alone it is 1.6 x 10 10 L & . The average 8 /im luminosity density at z < 0.3 
is 3.1 x 10 7 L & Mpc -3 , and the average luminosity density from dust alone is 
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2.5 x 10 7 Lq Mpc 3 . This luminosity arises predominantly from galaxies with 
8 /im luminosities (vL u ) between 2 x 10 9 and 2 x 10 10 L & , i.e., normal galaxies, 
not LIRGs or ULIRGs. 

Subject headings: cosmology: observations — galaxies:dust emission — galaxies: 
Survey — galaxies: mid- infrared 



1. INTRODUCTION 



Galaxy emission in the mid-infr ared (3/im-lO/im) arises from three main mechanisms 
(e.g. iPiovan. Tantalo. fc Chiosil 120061 ): 1) the Rayleigh - Jeans tail of stellar emission; 2) emis 



sion from heated dust i 


n the interstellar medium 


(ISM. Cesarskv et al. 1996 


Roelfsema et al. 


1996: 


Verstraete et al. 


1996: 


Vermeii et al. 


2002 


); and 3) a power law component powered 


by accreting black holes f 


Moorwood 


1986; 


Roche et al. 


1991; 


Elvis et al. 


1994 


). The ISM 



component displays prominent, broad emission features at 3.3, 6.2, 7.7, 8.6, and 11.3 /mi 
jGillett et all Il973l : Iwillner et all Il977l : IPhillips. Aitken. fe Rochel Il984h often attri buted 



1984 



Allamandola et al. 1985 



to polycyclic aromat ic hydrocarbons (PAH, iLeger fc Puget 

Puget fc Legerl 119891 ). Elliptical galaxies typically have little or no dust, and therefore 



their mid-infrared SEDs are expected to be dominated by stellar emission (ILu et al.l 12000 



20031 : IPahre et alil2004al lbl) . Star-forming galaxies, including the extreme Luminous Infrared 



10 12 Lm), show the PAH emission features ( 


Lu et al.ll2003 


:|Pahre et al.l2004a: 


Engelbracht et al. 


2005; 


Genzel et al. 


1998; 


Rieopoulou et al. 


1999; 


Armus et al. 


2007 


) . Active nuclei can occur 



in any galaxy type. Galaxy spectral energy distributions (SEDs) are the sum of radiation 
from all three mechanisms, and because the proportions vary widely, galaxy SEDs in the 
mid-infrared are very diverse. 

Mid-infrared dust emission luminosity is a good indicator of a galaxy's star formation 
rate. It is insensitive to dust obs curation, and it can b e measured for large samples of 
galaxies at both z ~ and z ~ 2. iRoche fc Aitkenl (119851 . and earlier papers) showed that 



mid-infrared feature emission is prominent in star forming galaxies, and ISO spectroscopic 



Vieroux et al. 


1999; 


Elbaz et al. 


2002). 



Forster Schrei 



jer et al. f|2004h 



Wu et all fl2005h 



rate (SFR — e.g. 

argued for a quantitative relationship between PAH strength and SFR, and 
confirmed a relationship using Spitzer IRAC 8 /im data to measure the PAH emission. 

Although PAH strength is a good SFR indicator, it is not perfect. Galaxies with 
powerful AGN may have very weak or no PAH emission features visible in their SEDs 
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( Moorwood 


1986 




Roche et al. 


1991 


2004 




Alonso-Herrero et ali 


2005 





19981 : lLaurent et all l2000l ; lEgami et al. 



This doesn't mean active star formation is absent be- 
cause, for examp le, intense ultraviolet ra diation produced by the accreting black hole may 
dissociate PAH (IRigopoulou et al.l Il999l ). especially near the nucleus. PAH features are 



also weak or absent in the SEDs of galaxies with low metallicities (lEngelbracht et al. 



2005 



Rosenberg. Ashby. Salzer. fc Huan32006 ; Wu et al.|[2006 ; Madden. Galliano. Jones, fc Sauvage 
20061 ) or low luminosities (IHogg et al.ll2005l ). At an extreme, iHouck et all (120041 ) found that 
the SED of the blue compact galaxy SBS 0335—052 shows no PAH emission at all. One 
possible explanation is that the formation and existence of PAH depends on metallicity. 
Whatever the limitations of PAH emission as an SFR indicator, it is still important to know 
how much PAH emission there is in the local Universe and what kinds of galaxies that emis- 
sion comes from. The results should at least give the SFR from typical galaxies, if not from 
unusual galaxy types. Understanding the local population is an essential step in determining 
how the SFR has evolved with time. 

The Spitzer Space Telescope offers the opportunity to observe large samples of galaxies 
at 8 /im. At low redshift, Spitzer's Infrared Array Camera (IRAC, Fazio et al. 2004) can 
detect the 6.2, 7.7, and 8.6 /im PAH emission features with its 5.8 and 8.0 /im detector arrays, 
and IRAC's wide field-of-view (5' x 5') and high sensitivity allow IRAC to map large areas 
of sky efficiently. Maps at 24 /im made with the Multi-band Imaging Photometer (MIPS - 



Rieke et al. 



Web 



200 4J) are an ideal w a y of selecting large samples of galax i es with PAH emission 



at z ~ 2 (IPapoyich et al 



d et al 



2005 



2006 



Lutz et al. 



Houck et al. 



2004 



Caputi et al 



200 



5; 



Houck et al 




2005 




Papovich e1 


: al. 


2006; 


authors (eg 


Houck et al. 


2005 




Yan et al. 



Lagache et al.l I2004J ) have used the Infrared Spectrograph (IRS - 
20041 ) to study the PAH emission features in the high-redshift 2 < z < 3 



population. 

This paper focuses on the 7.7 plus 8.6 /im PAH emission features from galaxies at low 
redshift detected in the IRAC 8.0 /im band. By combining redshifts and IRAC photometry, 
we obtain an IRAC 8.0 /im-selected local galaxy sample and derive a local PAH luminosity 
function. The IRAC 8.0 /im-selected galaxy sam ple is selected at th e same rest wavelength 



as a MIPS 24 /im-selected sample at z = 2 (e.g., ICaputi et al.l 120061 ) or an AKARI 16 /im- 
selected sample at z — 1. A comparison between these populations can thus reveal the effects 
of evolution between z = 2 and the present. The paper is arranged as follows: §2 describes 
the IRAC data and the sample selection. §3 discusses K-corrections for the IRAC bands 
and shows the color-luminosity relations of the sample galaxies. The luminosity function is 
presented in §4, and §5 summarizes the results. Distances throughout this paper are based 
on H = 70 km s" 1 Mpc" 1 , Q M = 0.3, Q A = 0.7. 
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OBSERVATIONS AND SAMPLE SELECTION 



The IRAC data come from the IRAC Shallow Survey (jEisenhardt et al.ll2004l ). which is 
one of the major IRAC GTO programs. The data cover 8.5 square degre es in the Bootes field 
of the NOAO Deep- Wide Field Survey (NDWFS Ijannuzi fc Devlll999l ). Each sky location 
was observed with three 30 s frames at 3.6, 4.5, 5.8, and 8.0 fim. To create the images in each 
waveband, the Bas ic Calibrated Data (BCD ) were mosaiced using the Spitzer Science Center 



(IBertin fc Arnouts 



software MOPEX (Makovoz fc Khanll2005l ). An 8 fim catalog was created using SExtractor 



19961 ). The 5a limit ing magnitude at 8.0 jum is 14.9 in the Vega system 



equivalent to a flux density of 69 /xJy (IBrodwin et a 
this limit is about 10600 galaxies deg -2 ( jFazio et al 



2004bf ). 



2007). The galaxy surface density at 



Photometry in all four IRAC bands was performed using the SExtracto r double-image 



mode . For each source in the 8 /im catalog, two magnitudes were derived (IBrodwin et al. 



20071 ): (1) an aperture magnitude within a 5" diameter, corrected to total magnitude using 
the point spread function (PSF) growth curve for point sources, and (2) a SExtractor "auto" 
magnitude. Figure [1] shows the difference between these two magnitudes as a function of the 
aperture magnitude. For objects with angular diameters larger than a few arcseconds, the 
aperture correction does not fully account for flux outside the aperture, and the aperture 
magnitude therefore underestimates the source flux. For faint point sources and also for 
extended sources with low surface brightness, SExtractor chooses too small a size for its 



"auto" aperture, and the auto magnitudes underestimate the flux (e.g., Fig. 5 of lLabbe et al. 



20051 ; iBrown et al.ll2007l ). In order to account for these biases, we used the brighter of the 
two magnitude measurements. For the largest galaxies the auto aperture encircles almost all 
of the flux and therefore gives total magnitude directly, while for faint galaxies, which are 
nearly pointlike, the aperture magnitude calibrated for point sources likewise lead to a total 
magnitude. Flux densities for low surface brightness galaxies may still be underestimated, 
but these are a small fraction of the sample^] All magnitudes in this paper are reported in 
the Vega magnitude system. 

In addition to th e photometric data, the AGN and Galaxy Evolution Survey (AGES, 
Kochanek et al.l 120061 ) has measur ed 15,052 redshifts i n the Bootes field with the Hec- 



Fabricant et al 



tospec multi-object s pectrograph (IFabricant et al.l Il998l ; iRoll. Fabricant. &: McLeodl [1998 



20051 ) on the MMT. Redshifts have been measured for 92% of galaxies to a 
magnitude limit of [8.0] = 13.2 and 65% of galaxies to [8.0] = 13.8 in a 6.88 deg 2 sub-field 
indicated in Figure [2j We limited the sample to this subfield and to [8.0] < 13.5 (251 /xJy) 



1 Less than 5% of galaxies have surface brightness lower than the equivalent of 13.5 mag in a 5" diameter. 
This is a strong upper limit on the fraction of galaxies that could be affected. 
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to minimize the redshift incompleteness and to z < 0.3 to minimize uncertainties in the 
K- correct ions E| Beyond z = 0.3, the strongest PAH feature at 7.7 /zm shifts out of the IRAC 
8.0 /im band. The sample was also limited to z > 0.02 to avoid statistical noise from small 
numbers of nearby galaxies. Figure [3] shows the redshift completeness in the final sample. 
The sample limit corresponds to a limiting 8 fim absolute magnitude of —27.7 at z = 0.3 
(neglecting the K-correction). 

The 8.0 /zm-selected sample includes normal galaxies, AGN, and stars. Figure H] shows 
the [5.8] — [8.0] color as a function of [3.6] — [4.5] for all of the objects in the sample. At 
low redshift, AGN and galaxies with PAH emission b oth have red IRAC colors, but th e 
two classes can be separated using IRAC color criteria (ILacy et al.l 12004 ; IStern et al.l 120051 ). 
Galactic stars and elliptical galaxies are mixed together at colors near (0,0); morphology was 
needed to distinguish between them. This sample is nearby enough that elliptical galaxies 
are detectably extended in the NDWFS data. The number of objects of each type are listed 
in Table 1 — 964 objects, 14% of the photometric sample, are identified as agnE However, 
only 42 of the 549 AGN with a measured redshift are at z < 0.3 whereas 2556 out of the 
3667 normal galaxies are at z < 0.3. The AGN fraction at z < 0.3 is therefore about 1.6% 
before allowing for incompleteness. Because AGES preferentially observed AGN candidates, 
the true AGN fraction is likely to be even lower than this. Known AGN were excluded from 
the sample, but separate results for them are explicitly mentioned in a few places below. As 
the results show, unknown AGN should have an insignificant effect on the derived luminosity 
functions. 



3. K-CORRECTIONS 

K-corrections for the galaxies in this sample are based on a simple two-component model 
for the SEDs. In this model, the SED for every galaxy is a linear combination of two fixed 
SEDs: (1) an old "early- type" stellar population, as might be found in an elliptical galaxy 
or spiral bulge, and (2) a mix of stars and interstellar emission as might be found in a "late- 
type" spiral galaxy disk. For each of these components, the SED was determined from an 
average of the observed SEDs for nearby galaxies that have both 2MASS photometry and 



2 The magnitude limits for the AGES spectroscopy, 15 < I < 20, had no significant effect on the sample 
selection. 

3 Some galaxies identified as "normal" probably have weak or obscured AGN, but these AGN are not 
detectable in either the IRAC colors or the AGES spectra and are therefore contribute little if anything to 
the observed light. 
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ISO spectroscopy (ILu et al.ll2003l ). Ellipticals were used to define the early-type SED, and 
disk galaxies were used to define the late- type SED0 The SED for each of these components 
was normalized at the rest-frame K band. Eleven model SEDs, ranging from a pure early- 
type template to a pure late- type template, were created. The late/early ratio for these 
model SEDs is defined as the ratio of the flux densities of the late-type and early-type 
components in the rest-frame K band. In addition, a twelfth template matching the M82 
SED was included to represent galaxies with high star-formation rates! Figure [5] shows the 
templates from which the model SEDs were constructed. Convolving each model with the 
IRAC filter functions gave the K-correction for each of the four IRAC bands as a function of 
redshift. Figure [6] shows the results. For z < 0.3, the 3.6 and 4.5 /im K-corrections are almost 
independent of the models, and the 8.0 /zm correction is also nearly model-independent unless 
the galaxy is early-type-dominated. In contrast, the 5.8 /im K-correction is highly variable 
because it depends critically on the strength of the 6.2 /im PAH feature. 

The suite of SED models spans the observed range of galaxy colors with very few outliers 
and almost none at z < 0.3. Figure [TJ shows the observed [4.5] — [8.0] galaxy colors as a 
function of redshift. Given the observed color and redshift for each galaxy, the proportion 
of early- and late-type SED components was determined by interpolation. This proportion 
defines the rest-frame colors for the galaxy and thus the K-correction for each of the IRAC 
magnitudes. The few galaxies with [4.5] — [8.0] colors redder than M82 were assigned the 
M82 K-corrections. AGN (when considered at all) were assigned K-corrections based on a 
power law SED between 5.8 and 8.0 jum. 



RESULTS 



4.1. Mid-Infrared Colors of the Sample Galaxies 



The galaxy sample divides into two populations in a color-magnitude diagram. Figure [7] 
shows the 8.0 jum absolute magnitude as a function of rest-frame [4.5] — [8.0] color. Galaxies 
in the red, late-type population ([4.5] — [8.0] ros t > 1.1, cf. Fig. [ 8|) are st ar-forming systems with 
prominent 8 /im PAH emission. According to iPahre et al.l (l2004bl ). this color corresponds 



4 We use "early-type" and "late-type" as convenient terms to refer to the two SED components despite 
having no direct knowledge of the actual morphology of individual galaxies. For example, irregular galaxies 
may have a "late-type SED" but not disk morphology. Mid-infrared colors for the SED components are the 
opposite of those in visible light: late types are red, and early types are blue. 



5 M82 has an IR luminosity of 4 x 10 10 L Q (|Belll 120031 ) corresponding to a SFR of about 7 M Q yr" 
SED is similar to SEDs of galaxies of even higher luminosity and SFR. 



Rs 
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to Sa and later morphological types. The blue, early-type population consists of galaxie s 



with weak or no PAH emission corresponding to SO/a and earlier types (IPahre et al.ll2004bl ). 
Figure M shows a different color-magnitude diagram, this one based on 3.6 /im magnitude, 
which is a measure of the stellar luminosity. The early-type population generally has higher 
stellar luminosity than the late-type population, and galaxies fainter than [3.6] a b s = —24 
are uncommon in the early- type population but common in the late- type one. Within the 
late-type population, the lack of correlation between the rest-frame [4.5] — [8.0] color and 
[3.6] a bs implies that the 3.6 fim luminosity is little affected by the 3.3 /im PAH emission 
feature. In contrast, [8.0] a bs is strongly correlated with the rest-frame [3.6] — [4.5] color 
(Fig. [10]), especially for [8.0] a b s < —27. This correlation suggests that dust emission at 
4.5 /im is detectable in the star-forming late-type population, i.e., galaxies with strong PAH 
emission also show dus t conti nuum emission at 4.5 /im. This is consistent with ISOPHOT 



spectroscopy: iLu et al.l (120031 ) detected 4 /im dust emission from nearby spiral galaxies and 
suggested it arises from a dust component closely related to the PAH. The same UV sources 
that heat PAH molecules may also heat the dust that produces the excess 4.5 /im emission. 
Dust emission at 4.5 /im is relatively unimportant for the least luminous galaxies, those 
with [8.0] a bs > —26, despite the presence of 8 /im emission in a subset of the low-luminosity 
galaxies (Fig. 0). 

The 8 /im number counts are strongly dominated by the late-type population: only 
3% of the spectroscopic sample have early-type colors. Even for the early-type-dominated 
galaxies, most of the 8 /im emission comes from dust rather than stellar light (Fig. [HD- The 
late-type galaxies are even more dominant in 8 /im luminosity, as shown in Fig. [7J All the 
early-type galaxies have 8 /im absolute magnitudes at least three magnitudes fainter than 
the brightest late-type galaxies. 



4.2. The 8.0 /im Luminosity Function 



There are many ways to estimate the galaxy luminosity function from survey data. 
When examining a new wavelength regime, a non-parametric method — one not assuming a 
specific functional form — is desirable because the true luminosity function may not follow a 
Schechter or other typical distribution. For the present survey, where galaxy distances vary 
by a facto r of ten, it is a lso b etter to avoid methods that are sen sitive to possible galaxy 
clustering. IWillmerl (119971 ) and lTakeuchi. Yoshikawa. fc Ishiil (120001 ) have presented detailed 
comparisons of many metho ds. We have compute d the lumino s ity fu nctions with both 
the familiar 1/V ma , x method JSchmidtlfl96~8l ) and the ILynden-Belll (ll97ll ) C method. The 
C~ method does not require binning the data, thus avoiding the possible bias in 1/Vmax found 
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by lPage fc Carreral (120001 ). and it is insensitive to possible galaxy clustering. The results of 
the two methods are consistent. A further verification that large scale structure has little 
effect was to ignore data in narrow but heavily populated redshift ranges such as z ~ 0.13 
and z ~ 0.19 (Fig. [HD- The resulting luminosity functions do not change significantly. 

The luminosity functions have to be corrected for spectroscopic incompleteness in the 
sample. Because this sample has spectra for nearly all galaxies (Table 1, Fig. [3]), the incom- 
pleteness correction has an insignificant effect on the results except in the faintest luminosity 
binsj^l This was checked by recomputing the luminosity function with two different incom- 
pleteness corrections. One method was to use only galaxies with spectroscopic redshifts but 
apply a weight to each one according to the reciprocal probability it would have been ob- 
served in the spectroscopic sample. The second method included all galaxies. Galaxies were 
divided into bins based on [4.5] — [8.0] colors and [8.0] magnitudes, and galaxies without 
redshifts were assigned the same redshift distribution as the galaxies with known redshifts in 
the same bin0 The two methods give luminosity functions that are equal within the Poisson 
uncertainties. 

The non-parametric 8.0 /im luminosity functions for various galaxy classes are shown 
in Figure [11] and Tables 2 and 3. As shown in this figure and also Figure 0, blue galaxies 
contribute little to the luminosity function and are rare for [8.0] a b s < —26. AGN are insignif- 
icant except in the brightest luminosity bins, where the uncertainties are large because few 
objects are in the survey. The luminosity function in the range —29.5 < [8.0] a b s < —22 is 
dominated by galaxies with significant dust contribution to the 8 /im flux. 

In addition to total luminosity, the observations allow us to measure the 8.0 /im PAH 
(dust) luminosity alone. For each galaxy in the sample, the rest-frame 3.6 Atm flux density 



was sc aled by a factor of 0.227 (an estimate of the stellar flux density at 8 /im — iPahre et al. 



2004al ) and subtracted from the rest-frame 8.0 /im flux density. This subtracts the starlight, 
leaving the PAH flux density. Because the subtraction is done on rest-frame flux densities, 
the procedure is equivalent to the two-component decomposition illustrated in Fig. [81 After 
subtracting the continuum, galaxies with residual 8.0 /im magnitudes fainter than the limit- 
ing magnitude [8.0] = 13.5 were excluded from the sample, and the luminosity function was 
re-calculated for the remaining galaxies to give a PAH-only luminosity function. The result 



incompleteness for bright apparent magnitudes — mostly galaxies at small distances — is distributed 
across all luminosity bins and has little effect on most of them because it affects only small numbers of 
galaxies in each bin. 

7 In practice, this was achieved by weighting the galaxies with known redshifts according to (total num- 
ber)/(number with known redshift) where the ratio was calculated for each bin. 
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is shown in Figure [TT] and Tables 2 and 3. 



The uncertainties in the luminosity functions are set both by Poisson statistics of the 
sample and by cosmic variance. The Poisson uncertainties for the l/V mauX method are given 
in Table 2 and statistical uncertainties for the C~ method in Table 3. Cosmic variance 



uncertainty can be estimated (IDavis fc Huchral [1982J) from the volume sampled plus an 
estimate of the galaxy two-point correlation function. This is not known directly for the 
8 jum sample, but we assume it's the same as for optically-selected galaxy samples. (If 
anything, the correlation function is likely to be smaller for 8 /im-selected galaxies, which 
ar e predominantly late - type, than for optically-selected ones, more of which are early-type 



Norberg et al. 



2002 



This will lead to less cosmic variance for an 8 /xm-selected sample.) 
In practice, a correl ation function based on the fluctuation power spe ctrum from WMAP-1 
( Spergel et al.ll2003h was extrapolated to the present via the method of lSeljak fc Zaldarriaga 
(119961 ) then transformed vi a a spherical Bessel func tion and smoothed with a 1 Mpc radius. 
The volume integral (e.g., iNewman fc David 120021 . Eq. 1) was then evaluated via a Monte 
Carlo approach. The resulting cosmic variance uncertainty is ~15% in bins [8.0] a b s < —27.7, 
where the survey samples the full volume of z < 0.3. This is a single uncertainty for 
the entire luminosity function, not an independent uncertainty in each bin. In the fainter 
bins, the uncertainty should in principle be larger, but the entire method is based on the 
assumption that the luminosity function is constant within the volume surveyed. In the 
faintest bin, for example, the survey magnitude limit corresponds to 135 Mpc, and the data 
give no information whatsoever on the luminosity function outside this distance. Density 
peaks in the sample are evident in Figure [H] but do not affect the derived luminosity function. 

As it turns out, a Schechter function is an excellent fit to both the total 8.0 fim and the 
PAH luminosity functions, as shown in Figure [TTJ Table 4 gives the best-fit paramet ers as 
found using the STY maximum likelihood method (jSandage. Tammann. fc Yahil | l979l ). and 
Figure [T21 shows the interdependence of the derived characteristic mag nitudes M% and faint- 
end slopes a. The derived Schechter exponents for the middle to faint end slopes are a = 
—1.38 and a = — 1.26 respectiv ely. The slopes are simila r to those found in the Sloan r band 



(ILovedav 



2003 



2004 



Xia et al 



2006 ) and the infrared K band (IHuang. Glazebrook. Cowie. &: Tinnev 



Jones. Peterson. Colless. fc Saundersl2006 ) over a similar redshift range, though Blanton et al 



(120031 ) found a slightly shallower slope in r. The steeper slope for the total luminosity shows 
that fainter galaxies tend to have relative less PAH emission compared to their stellar emis- 
sion (though the slope uncertainties overlap at the 2<r level, as seen in Fig. [12] and Table 4). 
The inactive "blue" galaxies have just as steep a faint-end slope (a = —1.39) as the entire 



8 The characteristic absolute magnitude values refer to luminosity vL v emerging from the galaxy near 
8 \xm measured in units of the Sun's bolometric luminosity. 
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sample, in contr ast to the result in blue light, where inactive g alaxies show a much shallower 
faint-end slope (IMadgwick. Hewett. Mortlock. fc Lahavl 120021 ) . At the bright end, L > L*, 
there are no e xcess galaxy numbers above the Schechter function prediction, in contrast to 
the results of ISaunders et al.l (119901 ) for the IRAS 60 /im luminosity function. The high 
luminosity populatio n, i.e., LIRGs and UL IRGs, consists of galaxies with far-infrared color 
temperatures >36 K (ISaunders et al.lll990l ). characteristic of Seyfert and starburst galaxies. 
Our exclusion of AGN in the 8.0 /im sample, which removed galaxies from the bright end of 
the luminosity function, does not explain why a Schechter function can fit the remainder of 
the sample. 

Integrating the luminosity functions over the whole sample gives the total 8 /im luminos- 
ity density in the local Universe. The result is 3.1 x 10 7 L & Mpc" 3 for total 8 /mi luminosity 
density (not including AGN), 2.5 x 10 7 L Q Mpc -3 for PAH luminosity density, and a highly 
uncertain 0.05 x 10 7 L & Mpc -3 for AGN luminosity density. (The sample is too small to 
determine reliable AGN numbers.) Thus over the entire galaxy population, stellar emission 
contributes about 19% of the 8.0 /im luminosity while dust emission contributes 80% and 
AGN roughly 1%. Galaxies with 2 x 10 9 L Q < uL v (8.0fjm) < 2 x 10 10 L contribute 50% 
of the total 8 /im luminosity, and galaxies above and below this range contribute about 25% 
each. Thus the dominant contribution to the PAH luminosity density comes from normal 
galaxies, not LIRGs or ULIRGs. 



4.3. Star Formation Rates 



Star formatio n rates in the local universe have be en measured by a variety of techniques. 
(See the review by lKennicuttlll998al ) iHopkinsi (120041 ) compiled star formation rate densities 
(SFRD) from a variety of tracers. According to Hopkins' compilation, over the z < 0.3 
range of the present survey, the volume-weighted average SFRD « 0.031 M & yr _1 Mpc~ 3 , 
though some tr acers tend to be syst ematically higher (emission lines) or lower (ultraviolet) 
than this value. iMartin et al.l (120051 ) combined ultraviolet with far infrared data and found 
a local (z < 0.04) SFRD « 0.021 M yr" 1 Mpc" 3 , about 20% higher than but consistent 
with Hopkins' result at this much smaller redshift. 

Dust emission at 8.0 /im is correlated with the total infrared luminosity Lm in galaxies 



(Rig opou 



(ILu et al 



ou et al. 1999 ; Elbaz et al. 2002 ; Lu et al. 20031) although the rel ation is non-linear 



2003) and has large scatter (e.g., Fig. 17 of ISmith et al.l 120071 ). In normal star- 



forming g alaxies, the integrated flux betw een 5.8 an d 11.3 /im is between 9% (fo r galaxies 



Lu et al.l call "FIR-quiescent" ) and 5% (for lLu et al.l "FIR-active" g alaxies) of Ltr (ILu et al. 
20031 ) . Recent results from Spitzer are consistent with these values; Smith et al. ( 2007 ) find 
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that in the central areas of normal galaxies, the 7.7 fim PAH complex can contribute up to 
10% of Ljr. In ULIRGs, PAH emission in the 5.8-11.3 /im range is lower at <1% of Ljr 
( IRigopoulou et al.lll999l ; lArmus et al.ll2007l ) ; this coul d be partly due to an AGN contributing 



to Lj R . Despite the uncertainties, some authors (e.g. jForster Schreiber et al.ll2004l ; IWu et al. 



20051) have suggested that vL u (PXH) can be used to measure star formation rates for galaxies, 
and IWu et ahf fl2005h derived SFR « i/L 1/ (PAH)/(1.5 x 10 9 L ) M yr" 1 . This result was 
based on a sample of Sloan galaxies show ing Ha emissio n. The PAH luminosity density found 
in §4.21 combined with the local SFRD (jHopkins! 12004 . Fig. 2) integrated from z = to 0.3 
gives a volume average SFR « vL v (PAH)/ (0.84 x 10 9 L ) M Q yr" 1 , about double that found 
by IWu et al.l for individual galaxies. In other words, the volume- averaged PAH luminosity 
density is half as much as expected, given the volume-averaged SFRD and SFRD/i/L^(PAH) 
ratios of individual galaxies. The difference suggests that the galaxy samples have been 
weighted to galaxies with unusually large PAH emission for their SFRs, but this is hard to 
understand if th e dominant contribut ion to star formation comes from the normal galaxy 
population (e.g.. lLe Floc'h et al.ll2005l ). 



5. SUMMARY 

An 8.0 /xm-selected sample of low redshift (z < 0.3) galaxies finds primarily star-forming 
galaxies. AGN represent less than 2% of the sample both by number and by fraction of the 
8 (Ltm luminosity. The luminosity function of the remaining galaxies can be described as a 
Schechter function with characteristic magnitude M* = —28.46 (corresponding to vL v = 
1.8 x 10 10 Lq) and a = —1.38. Subtracting starlight continuum emission from each galaxy 
gives a luminosity function for emission from the interstellar medium component (mostly 
PAH) alone. This can also be approximated by a Schechter function with characteristic 
magnitude M* = -28.34 (uL u = 1.6 x 10 10 L & ) and a = -1.26. 

The total 8.0 /im luminosity density for this 8.0 /im-selected sample is 3.1 x 10 7 L & Mpc -3 . 
Ignoring the luminosity contributed by AGN, about 81% of the 8 /im luminosity is from in- 
terstellar dust emission, presumably PAH, and about 19% is from stellar photospheres. The 
observed, volume-averaged ratio of PAH luminosity density to star formation rate density is 
about half of an estimate based on local, individual galaxies. 
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by JPL / Caltech. JLR has received support from an NSF Astronomy and Astrophysics Post- 
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Fig. 1. — Difference between the SExtracter total ("auto") magnitudes and 5" aperture 
magnitudes. Offsets are mostly positive, indicating that aperture magnitudes underestimate 
the total flux of resolved galaxies. We selected the sample using the aperture magnitudes to 
avoid bias against lower surface brightness galaxies, and set the magnitude limit at [8.0] = 
13.5. The derived luminosity functions are based on the auto magnitudes for the resolved 
galaxies and the aperture magnitudes for the unresolved galaxies. 
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Fig. 2. — The Bootes field showing the approximate area of 8 //m coverage outlined in red. 
Positions are J2000. Black dots denote denote the 8.0 /xm-detected galaxies that have AGES 
spectroscopic redshifts. The black lines enclose the area of the highly complete redshift 
sample used for this paper. 
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Fig. 3. — Redshift identification fraction in each half-magnitude bin (filled circles, left ordi- 
nate label). Error bars show the Poisson uncertainty. The solid line shows the histogram of 
the 8 /im magnitude distribution of all galaxies in the catalog (right ordinate label). Only 
galaxies with [8.0] < 13.5 are in the final sample. 
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Fig. 4. — [5.8] — [8.0] vs [3.6] — [4 .5] color-color dia gram for separation of AGN, star-forming 
galaxies, and early type galaxies (IStern et al.ll2005l ). Data are plotted as observed, i.e., before 
K-correction. Black dots denote normal galaxies, red dots denote stars, and blue dots denote 
AGN, all as classified by optical spectroscopy. The two thick lines show two SED models 
from < z < 0.6, one with and the other without PAH emission. (T he latter is almost 
hidden behind the data near zero color.) The thin solid lines show the IStern et al.l (120051 ) 
AGN-galaxy separation criterion; AGN lie within the quadrilateral extending to the top of 
the plot. The color cut at [5.8] — [8.0] = 0.5 separates stars and elliptical galaxies from 
galaxies with significant PAH emission. Galaxies with PAH emission cannot enter the "no 
PAH" area unless their redshifts are >0.6. 
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Fig. 5. — Galaxy SED templates. The upper curve shows the M82 template, the middle 
curve the spiral ( "late- type" ) template, and the bottom curve the elliptical ( "early- type" ) 
template. The vertical axis is in arbitrary units. The early-type te mplate is an ave rage 
ISOPHOT-S spectrum of six classical elliptical gal axies compiled by iLu et al.l (120001 ) . It 
is very similar to the elliptical spectrum shown by ILu et al.l (120031 ) except at wavelengths 
longer than 8 /im, where an improved ISOPHOT-S calibration made the current spectrum 
slightly s t eeper (i.e., bluer). The late-type template is the average for spiral galaxies from 
Lu et al.l (120001 ) . Nine additional templates used in this paper are mixtures of the early- and 
late-type templates. 



-22 - 



1 .0 




Redshift Redshift 

Fig. 6. — K-corrections for the four IRAC channels for each of the twelve SED models. 
The bottom line in each panel shows the early-type model. The 3.6 and 4.5 /im bands 
sample mainly stellar emission longward of the Wien peak and therefore show a negative 
K-correction. The 5.8 /im band is relatively narrow and shows a positive K-correction, 
depending on the amount of PAH present, as the 6.2 /im PAH feature leaves the band at 
z ~ 0.05. As redshift increases further, the band samples only stellar emission, and the 
K-correction decreases. The 8 /im band is wide, and the 7.7 /im feature remains inside 
the band at low redshift. As redshift increases, the 6.2 /im PAH feature enters the band 
before the 7.7 /im feature leaves, leading to a nearly constant K-correction. At z > 0.2, the 
7.7 /im feature begins to leave the band, and the K-correction increases according to the 
PAH strength. 
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Fig. 7. — Color-magnitude diagram showing [4.5] — [8.0] color versus 8 /im absolute magni- 
tudes. Points denote galaxies with spectroscopic redshifts. Colors and magnitudes are both 
corrected to the rest frame. Most galaxies show strong PAH emission ([4.5] — [8.0] > 2) 
and are therefore in the "late- type" category as we have defined it. "Early- type" galaxies 
with weak or no PAH emission ([4.5] — [8.0] < 0.4) are well separated in this diagram from 
the much larger late-type population, but galaxies span the entire color range. None of the 
early-type galaxies has high luminosity at 8 jum. 
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Fig. 8. — The [4.5] — [8.0] color versus redshift for galaxies with spectroscopic redshifts. Points denote 
the observed galaxies with AGN excluded. Lines show colors for each of twelve model SEDs as a function 
of redshift. The reddest model is derived from the M 82 SED. Other models are a linear combination of 
late-type and early-type templates as shown in Fig. [5l These SEDs are normalized at rest 2.2 /im with the 
stellar fraction going from to 1 in steps of 0.1 from the bluest model to the reddest one (second-reddest 
including the M82 model). Only galaxies with z < 0.3 are used in this paper. 
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Fig. 9. — Color- magnitude diagram showing [4.5] — [8.0] color vs 3.6 /im absolute mag- 
nitudes. Points denote galaxies with spectroscopic redshifts. Colors and magnitudes are 
both corrected to the rest frame. The upper luminosity envelope is populated by both blue 
([4.5] — [8.0] < 1.1) and star-forming galaxies, but there are no blue galaxies with low 3.6 /im 
luminosities. For the strongly star-forming galaxies, [4.5] — [8.0] > 2, color does not correlate 
with the absolute 3.6 jum magnitude. 
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Fig. 10. — Color-magnitude diagram showing [3.6] — [4.5] color vs 8 /im absolute mag- 
nitudes. Points denote galaxies with spectroscopic redshifts. Colors and magnitudes are 
both corrected to the rest frame. The correlation between color and magnitude indicates 
that dust emission may constitute about 30% of the total 4.5/zm emission for galaxies with 
[8.0] abs < -27. 
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Fig. 11. — Derived differential luminosity functions corrected for spectroscopic incompleteness. Black 
dots denote the 8 /itm luminosity function for all galaxies excluding AGN, and red squares denote the 
luminosity function for PAH alone, i.e., after subtracting the stellar continuum emission for each galaxy. 
These two luminosity functions were derived from the non-parametric C - method. Blue triangles denote the 
luminosity function for early-type galaxies, specifically those with [4.5] — [8.0] < 1.1. Even for these galaxies, 
much of the emission is PAH despite the color selection. Green crosses denote the luminosity function for 
AGN. These latter two luminosity functions were derived from the 1/V max method because of the smaller 
numbers of galaxies. Because of the small numbers, the blue galaxy luminosity function is only reliable for 
—27 < [8.0] a h s < —25. The AGN luminosity function is highly uncertain throughout; only one bin has as 
many as 8 galaxies in it. All points are slightly offset horizontally for clarity; see Tables 2 and 3 for actual 
magnitude bins. Error bars show the Poisson uncertainties in each bin. Lines show the respective parametric 
Schechter functions that best fit the data for "all 8 /im emission" and "PAH only." Both were derived from 
the STY method. The bottom axis is labeled in 8 /im Vega magnitudes, while the top axis is labeled in the 
equivalent vL v units. 
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Fig. 12. — Uncertainty ellipses for characteristic magnitude M* and faint end slope a of 
the Schechter luminosity functions. The ellipses were calculated using a grid in M*-a space 
and the STY formalism. The curves represent contours of constant A% 2 corresponding to 
la (inner), 2a (middle), and 3a (outer) confidence levels assuming three fitting parameters. 
Statistical uncertainty in the overall normalization <£>* is very small, but the true uncertainty 
is set by cosmic variance (§4.2). Faint end slopes probably differ but are consistent at the 
2a level. 
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Table 1 

Number of Objects at [8.0] < 13.5 



Object Total Redshift Sample z < 0.3 

Star 1894 

Blue Galaxies 236 109 100 

Star-Forming Galaxies 3667 3174 2556 

AGN 964 549 42 



-30- 



Table 2 

Luminosity Functions Calculated by the 1/V max Method 



[8.0] afes a uL u (8 fim) N 0(8.O//m) A0(8.O//m) N 0(PAH) A0(PAH) 













10~ 3 Mpc" 3 
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Absolute Vega magnitude bin centers. 
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Table 3 

Luminosity Functions Calculated by the C~ Method 



[8.0] a6s a uL u (8 urn) 0(8.O//m) A0(8.O//m) 0(PAH) A0(PAH) 
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a Approximate bin centers in absolute Vega magnitudes. The actual first bin center is —30.15 for the "8.0" 
columns and —30.11 for the "PAH" columns. Successive bin centers differ by half- magnitude intervals. 
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Table 4 

schechter function parameters 





all 8 /xm emission 


PAH only 


M*(8 /xm) 


-28.46 ±0.08 


-28.34 ±0.08 


L*(8 p) 


1.8 x 10 10 


1.6 x 10 10 L Q 


$* a 


1.10 x 10~ 3 


1.17 x 10" 3 


a 


-1.38 ±0.04 


-1.26 ±0.04 


Integrated luminosity b 


3.1 x 10 7 


2.5 x 10 7 L Q Mpc~ 3 



"Fiducial galaxy density in galaxies Mpc 3 mag 1 . Statistical uncertainty is very small, but the real uncer- 
tainty is about 15% set by cosmic variance. 

integral of 4>(L)dL, where L = vL v for v corresponding to 8 /um. 
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